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ABSTRACT 

The probability that an existing planetary transit is detectable in one's data is sensitively dependent 
upon the window function of the observations. We quantitatively characterize and provide visualiza- 
tions of the dependence of this probability as a function of orbital period upon several observing 
strategy and astrophysical parameters, such as length of observing run, observing cadence, length of 
night, transit duration and depth, and the minimum number of sampled transits. The ability to detect 
a transit is directly related to the intrinsic noise of the observations. In our simulations of observa- 
tional window functions, we explicitly address non-correlated (gaussian or white) noise and correlated 
(red) noise and discuss how these two noise components affect transit detectability in fundamentally 
different manners, especially for long periods and/or small transit depths. We furthermore discuss the 
consequence of competing effects on transit detectability, elaborate on measures of observing strate- 
gies, and examine the projected efficiency of different transit survey scenarios with respect to certain 
regions of parameter space. 

Subject headings: eclipses — methods: statistical — planetary systems — surveys — techniques: 
photometric — time 



1. INTRODUCTION 



The signal-to-noise ratio (SNR) of a planetary tran- 
sit detection in photometric time series data can, in the 
simplest case, be approximated by: 
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noise, due to their focus on bright targets and high num- 
ber of observational epochs: both are aimed to reduce 
white noise, and therefore make red noise the dominant 
component. 

A detailed description of the transit detection SNR 
which in cludes both white and red noise components is 
given bv lPont et all (|2006l ): 



In this equation, depth is the transit depth in magni- 
tudes, a represents the photometric measurement uncer- 
tainty in magnitudes per data point (assumed here to 
be the same for all data points), and n equals t he to- 
tal nu mber of data points observed during transit (|Pontj 
2006). One essential assumption in this equation is the 
absence of any statistically correlated (red) noise, i.e., 
only random (white) noise is present. White noise is de- 
fined as noise that is uncorrelated from data point to data 
point; typical sources are photon noise and sky back- 
ground noise. The relative contribution of white noise 
to the total noise decreases with increasing brightness of 
t he obs e rved t arget and numbe r of da ta points. 

iPontJ (120061 ) and iPont et alj (|2006f ) showed that cal- 
culations of transit SNRs with only white noise, as in 
Equation[TJ is often insufficient and overly optimistic. In- 
stead, one needs to account for the presence of red noise 
in calculations of SNRs and the corresponding yield pro- 
jections for transit surveys. Red noise is defined as noise 
that is correlated from data point to data point; it is 
not necessarily removed through standard differential or 
ensemble photometry techniques. Typical sources of red 
noise may be weather, seeing changes, tracking/guiding 
errors, flatfielding errors, changes in airmass, or intrinsic 
astrophysical changes in target brightness. It does not 
change as a function of target magnitude, and is gen- 
erally independent of the number of obse rvational data 
points (see equation 9 in lPont et al.l [20061 ). Thus, plan- 
etary transit searches are particularly sensitive to red 
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where cov[i;j] is the covariance matrix with elements 
Cij representing the correlation coefficents between the 
i-th and j-th measurements obtained during transit. All 
diagonal elements Cu = of are not correlated with 
other measurements and thus represent the uncorrelated 
or white noise uncertainties in the i-th. measurement. 
These diagonal elements are assumed to be the same, 
i.e., <Ji = a for all values of i. 

In order to make th e abov e equation more practically 
calculable. iPont et al.l (|2006f ) assume that statistical cor- 
relation among data points from different transits will be 
much weaker than among data points observed during 
the same transit. They furthermore separate the total 
noise into a purely uncorrelated (white) component a w 
and a purely correlated (red) component a r and derive 
an approximation of Equation O 



SNR 



transit — 
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where n is the total number of data points observed dur- 
ing all transits, N tr is the total number of transits ob- 
served, rik is the number of data points observed during 
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the k-th transit, and a w and oy are the white and red 
noise components, respectively. 

By means of EquationsQ]and[31 it is clear that a planet 
transit SNR can be regarded as a function of transit sur- 
vey strategy and astrophysical parameters (see $2]). If 
this SNR exceeds a certain threshold value, then an ex- 
isting transiting planet is, for the purposes of this paper, 
defined to be detectable in the data 1 . 

The window function determines the probability, as 
a function of planetary orbital period, that SNR tra „ S it 
exceeds this threshold. In this paper, we examine the 
dependence of the detection probability upon several as- 
trophysical and transit survey strategy parameters for 
a number of white noise and red noise assumptions as 
well as criteria based on a minimum number of transits 
sampled. 

Since our calculations are based on existing transits, 
we note that the following aspects are not taken into ac- 
count: the estimated frequency of transiting exoplanets, 
any non-circular orbits, multi planet systems, and detec- 
tion of secondary transits. We also do not address the 
problem of false positives and how to weed them out. For 
more detailed studies of the above, we refer the reader 
to the f ollowing stud ies: freque ncy of (transiting) exo - 
planets: iGould et all (|2006D and lCumming et all (|2008l ): 
transit proba bi lity as a fun ction of orbital elements 
iBarnes! (|2007f ). iBurkd (120081) . a nd iKane fe von Braun 



(20C 



2008); plus see iGaudil 12007) and iBeattv fe Gaudi 



(2008) and references therein for a comprehensive study 
of all factors influencing planet detections in transit sur- 
veys. 

We briefly outline our methods in $2] which describes 
our algorithm in ^2.1| along with a justification for the 
threshold SNR selection, and addresses the respective 
influences of varying white and red noise components 
( §2.2p . as well as the consideration of sampling at least 
Ntr transits with one's data to constitute a detection 
f £|2.3[) . We examine the effects of various survey strategy 
and astrophysical parameters in Sj3j Section [4] contains 
the application of window functions for selected scenar- 
ios and types of survey. We summarize and conclude in 

m 

2. ALGORITHM AND PARAMETERS 

In this Section, we provide a brief description of our 
algorithm and explain our choices of the globally used 
values of input parameters in our calculations in |£3j 



2.1. Description of the Algorithm 

The window function algorithm used in this paper is 
based on counting data points observed during transit 
whose contribution to a virtual detection is dependent 
on the values of o~ w and oy as defined in Equations [2] 

1 We note that, to maximize applicability for astronomical 
planet transit sur v eys, w e follow the arguments outlined above 
and in Pont ct al. (2006), rather than using more rigorous treat- 
ments employed in the large body of statistical literature devoted 
to time-series analysis. These treatments include autoregressive 
moving-average (ARMA) modeling where events such as eclipses 
in the presence of red noise are sought with the help of autocor- 
relation functions at d i fferent time lag s l|Kocn & Lombard! 119931 : 
[Lombard & Kocn 1993; Robinson 2005|), power spectrum analysis 
(Konig & Thinner 199'if), or the use of surrogate data sets (JTimmerl 
119981) . 



and [3l typically measurable or calculable quantities in 
photometric time series surveys. 

User-provided observing cadence, number of nights, 
and typical length of night are used to generate an ob- 
serving time line. From the input stellar and planetary 
radii, we calculate transit depth and duration acco rding 
to the equations in Se ager fc Mallen-Ornelasl (|2003l ) (ex- 
cept in £13.61 where we explicitly set transit depth and du- 
ration), thereby assuming a central transit (i.e., impact 
parameter 6 = 0) and zero-length ingress and egress. For 
each orbital period, a family of light curves is generated 
for a range of starting phase angles; each with transits of 
user-defined photometric depths at the appropriate in- 
tervals. 

In the simulations, the number of data points per tran- 
sit (n&), number of transits (N tr ) and total number of 
data points within all transits (n) are tracked. It should 
be noted that an observation has to fit fully within a 
transit to be counted toward n and rik (that is, it needs 
to start after the beginning of the transit and terminate 
before the end of the transit), resulting in shorter expo- 
sure times' being more favorable for transit detection in 
this algorithm. For every light curve, the SNR (Equa- 
tions [T] and [3]) is calculated. If, for a given phase angle, 
the SNR exceeds SNR t hreshoid, a transit is considered 
"detected". The probability of detection (P detection) for 
a given orbital period is simply the ratio of phase angles 
for which a transit was detected to the total number of 
phase angles. 

Typical observational parameter values assumed in this 
paper (unless specifically noted) are: minutes for the ob- 
serving cadence, one minute for the exposure time, tens 
of nights for observing run length, and few to ten hours 
for the typical time of observation spent during one night 
on the monitored target. Astrophysical parameter values 
are assumed to be around 1.0 and 0.1 solar radii for the 
parent star and orbiting planet, respectively, resulting in 
a transit depth of 0.01 mag. Transit duration depends 
on period, but typical duty cycles are in the 1% to few % 
range. Additionally, we set o~ w and oy to a few millimag- 
nitudes (mmag). T he threshold SNR is set to 7.0, based 
on the arguments inlJenkins et al.l (|2002t ) and specifically 
iPont et al.l(l200ll2007h . which each use thresholds of 7- 
9 as acceptable values for reducing false-alarms whilst 
maximizing real detections given a typical transit survey 
configuration. 

We note that, in contrast to other some window func- 
tion calculations in the astronomical literature, we only 
use the SNR criterion to quantify detections, along with 
an assumed minimum number of sampled transits, and 
we do not require that, e.g. , a full transit be contained 
in the data (as in, e .g., iMallen-Ornelas et al.l 120031 : 
Ivon Braun et al.ll2005[ ). We do not account for holes in 
the observing due to weather, telescope outages, or tech- 
nical problems. Furthermore, as mentioned in SjTJ we 
only calculate the probability of detecting existiting pri- 
mary transits in circular orbits. Finally, we assume that 
the number of out-of-transit data points sampled is much 
higher than the number of in-transit data points. 

2.2. Red Noise and White Noise 



As ou tlined in IPepper fc Gaudil (120051): IPont et al" 
(2006); lAigrain fc PontJ (|2007f ): | von Braun fc Ciardil 



2008), red noise is the dominant source of noise in the 
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regimes of brightness and number of observational epochs 
where transit surveys typically operate since red noise is 
independ ent of target brig htness. We show this effect in 
SETT1 See lAigrain fc Pont) (|2007t ) and llrwin et~all (|2007f ) 
for an in-depth discussion of different noise properties 
and their calculations. 

Typical ground-based survey estimates of ay, as 
defined in £JTJ are on the order of 2-6 millimag- 
nitudcs (mm ag) (e.g iPont et all 120061: llrwin et all 
120071: iKane et alj|2008t iNutzman fc Charbonneaij 2008) . 
When subjected to detrending algorithms such as TFA 



TABLE 1 

Mean Pdetection Values for Different Magnitudes of o> 



Kovacs et al 



Tamuz et al 



[20051 : iKovacs fc Bakosll2"007h or SYSREM 
ay can be reduced to 1-2 mmag. 
It is worth pointing out that the influence of red 
noise is much less of a problem for targeted observations 
such as characteri zation of known planetary transits (see 

iGillon et al.ll2008l for e xample). 

Studies to date (e.g.. lAigrain et aLll2008f ) have shown 
that the red noise in the space-based CoRoT mission 
(Baglin et al. 2006) is significantly lower than in ground- 
based counterparts, due in large part to the "removal 
of the atmosphere" (|Pont et all 120061 : iBeattv fe Gaudil 
2008). Thus, the cause of any space based red noise not 
due to stellar variations is most likely caused by varia- 
tions in the thermal environment of the spacecraft and 
detectors. Typical values for o> in CoRoT light curves 
are on the order of 0.5 mmag (R. Alonso 2 008, private 
communication; see also lAigrain et aill2009f ). 

2.3. Number of Sampled Transits 

One criterion often used to calculate detection effi- 
ciency and related survey yield is the minimum number 
of sampled transits (i.e., the minimum number of tran- 
sits during which any data were obtained). An important 
factor in the succes s of the widely used BLS algorithm 
(jKovacs et al.1120021 ) is the initial folding of the data by 
a test period and subsequent search for transit-like fea- 
tures in the phased data. Thus, its power is really only 
realized for data that contain more than one sampled 
transit. We assume in this publication that the BLS al- 
gorithm has become an "industry standard" in the search 
for planetary transits, and we thus require the existence 
of at least two transits in the data for a transit detec- 
tion, except for where we explicitly change this criterion 
( £!3.5p . It is worth noting that different simulations in the 
literature require different minimum number s of transits 
sampled, such as three for lPont et al.l (|2006l ). 

3. THE INFLUENCE OF OBSERVATIONAL WINDOW 
FUNCTIONS AND ASTROPHYSICAL PARAMETERS ON 
TRANSIT DETECTION PROBABILITY 

Careful consideration of the various strategy as- 
pects involved in planetary transit surveys and a 
number of astrophysical parameters will have signif- 
icant eff ects on the detection effic i ency of existing: 
trans it s (iMallen-Ornelas et al.l |2003t Ivon Braun et al.l 
20051 IPepper fe Gaudil 120051: IBeattv fc Gaudil 120081: 
von Braun fe Ciardil 120081 : IBeattv! l2009ft . This Section 
quantitatively illustrates these effects under considera- 
tion of the assumptions described in ^\ For the sake of 
clarity, we vary one parameter at a time, leaving all oth- 
ers fixed to values justified in fjl] and <|2] In particular, 
.11 examines different values for red and white noise, 
looks at various observing run lengths (with a given 





oy = 


1 mmag 


2 mmag 


3 mmag 


4 mmag 


0-5 days 


0.983 


0.982 


0.978 


0.906 


0.540 


5-10 days 


0.945 


0.942 


0.819 


0.102 


0.010 


10-20 days 


0.753 


0.747 


0.235 


0.012 


0.000 


20-30 days 


0.425 


0.421 


0.051 


0.000 


0.000 



Note. — Mean values for Pdetection for various period ranges 
(column 1) as a function of different values of oy (Fig. [T}. As- 
sumed pa ram eters are given in the caption of Fig. [l] For discus- 
sion, see £13.11 

number of hours of observing per night), whereas £13.31 
assumes a number of consecutive observing nights but 
varies their lengths. £13.41 investigates different observing 
cadences. In £j33J we explicitly change the criterion of 
a minimum of two transits sampled for a detection that 
we mention in ^2.3l to see how requiring a larger number 
decreases detection efficiency. £13.61 deals with different 
transit depths and durations. 

The values of the parameters held constant in the re- 
spective calculation are given in the caption of the ap- 
propriate figure. The solid (blue) line shows the detec- 
tion efficiency in the hypothetical case of zero red noise, 
and the dashed (red) line shows the same for a given 
ay 7^ 0. The corresponding table shows the mean val- 
ues of Pdetection for the ranges of orbital periods given 
in the first column under assumptions of the indicated 
magnitudes of o~ w and ay. 

3.1. Amount of Red Noise and White Noise 

The contribution of red noise is independent of target 
brightness (unlike white noise, which is mostly due to 
photon noise for the brightest targets). Since planet tran- 
sits are typically detected around the brightest sources 
in a given dat a set, red noise will be the dominant 
source of noise (IPepper fc Gaudil 120051: IPont et alJ feOOe; 
lAigrain fc Pcmt] 120071: IBeattv fc Gaudil l2008f l. Figures [T] 
and [2] along with Tables [T] and [2] quantitatively substan- 
tiate this statement, illustrating the influences of differ- 
ent amounts of red and white noises for different period 
ranges. As we mention in ^2.21 typical values for the 
wide-field ground-based transit surveys before detrend- 
ing are o~ w ~ 5 mmag and ay ~ 2-6 mmag, which reduces 
to ay ~ 1-2 mmag after detrending. 

The difference in Pdetection in Fig. Q] and Table Q] be- 
tween ay = 1 mmag and ay = 4 mmag is very significant 
for longer periods. In addition, Fig. [2] and Table [2] show 
how small the influence of a w upon Pdetection is for no 
or very little red noise. Thus, the value of minimizing 
the influences of red noise during observing (even at the 
expense of increasing o~ w if necessary), an d of applying 
detrending al gorithms such as S YSREM dTamuz et al.l 
120051 ) or TFA (jKovacs et al.ll2005l : lKovacs fc Bakosl2007l) 
to one's data as part of their reduction can hardly be 
overstated. 

3.2. Observing Run Length 

For any kind of transit survey that has limited access 
to telescope time, the question of how long to spend on 
one field will occur at some point during the design of 
the observing strategy. At what point is it worth switch- 
ing to a different field to increase the number of targets 
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Fig. 1. — Variable red noise and constant white noise (a w = 5 mmag). The (blue) solid line indicates the detection efficiency for no red 
noise (a r = 0) and is shown in all four panels. The (red) dashed line indicates detection efficiency for varying levels of oy as indicated 
in the respective panel. Additional parameters are: Subthreshold = 7.0, 60 consecutive observing nights with 8 hours of observing every 
night, an observing cadence of 5 minutes, R s ta r = Rq, Rpianet = 0.1-Rq, M s tar = Mq, M p i anet <C Mq. See Table[T]for mean values of 
Pdetection over various period ranges, and £13.11 for details. 



TABLE 2 

MEAN Pdetection VALUES FOR DIFFERENT MAGNITUDES 
OF a w 





1 mmag 


3 mmag 


5 mmag 


10 mmag 




(o> = 0) 


0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.984 
0.947 
0.755 
0.428 


0.983 
0.947 
0.755 
0.428 


0.983 
0.945 
0.753 
0.425 


0.977 
0.916 
0.677 
0.355 






- = 2mmag) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.981 
0.897 
0.412 
0.131 


0.980 
0.857 
0.281 
0.058 


0.978 
0.819 
0.235 
0.051 


0.967 
0.444 
0.072 
0.017 



Note. — Mean values for Pdetection f° r various period 
ranges (column 1) as function of a w (Fig. |2j. Assumed 
parameters are g iven in the caption of Fig. [2] For dis- 
cussion, sec i|3.1| 



without overly reducing the probability of detecting ex- 
isting transits in the data? We provide insight into the 
answer to this question in Fig. [3] and Table EB 

To first order, observing a field for few nights will yield 
an almost negligible probability of detection, potentially 



leading to a waste of telescope time. Alternatively, it may 
be wise not to stay on a single field for too long but rather 
to double the chances of detecting any planetary transits 
by switching fields and thus increasing the number of 
monitored stars. It is ultimately a question of the period 
range one is sampling in a given survey. As Fig. [3] shows 
for a typical set of parameters, "very hot Jupiters", i.e., 
planets with perio ds up to ~3 days per the definition in 
iGaudi et all pTOl . can be detected even with a residual 
presence of red noise and "only" 15 nights (eight hours 
per night) of monitoring. However, longer period planets 
(~6 days and longer) remain elusive (for ay > 2mmag) 
until the length of the observing run exceeds 30 nights. 

3.3. Length of Night 

The amount of time for which a given target field can 
be observed from the ground during one night depends on 
its celestial coordinates, the location of the telescope, the 
time of year, and, of course, outages due to weather, or 
technical or other problems. Special cases are discussed 
below, such as space-based observing ( §4.1j) or synoptic 
surveys ( 34.20 . 

The length of night can also depend on observing strat- 
egy. As an alternative to decreasing the number of nights 
spent on a single target field to increase the number of 
monitored stars ( 33.2[) , one may instead choose to split 
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period[days] period[days] 



Fig. 2. — Variable white noise and constant red noise. The value for a w is given in the top right of every respective panel. The (blue) 
solid line indicates the detection efficiency a r = 0, and the (red) dashed for oy = 2 mmag. Additional parameters are: SNR^reshoid = 7.0, 
60 consecutive observing nights with 8 hours of observing every night, an observing cadence of 5 minutes, R s t ar = Rq, Rpianet = 0.1-Rq, 
M s tar = Mq, M p i ane t <S Mq. See Table [2] for mean values of Pdetection over various period ranges, and i|3,ll for discussion. 



TABLE 3 

MEAN Pdetection VALUES FOR DIFFERENT 

Observing Run Lengths 



Nights 


15 


30 


60 


90 






(cr, 


= 0) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.792 
0.238 
0.027 
0.000 


0.957 
0.696 
0.279 
0.064 


0.983 
0.945 
0.753 
0.425 


0.988 
0.969 
0.918 
0.736 




(oy = 


2mmag) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.482 
0.020 
0.000 
0.000 


0.865 
0.189 
0.022 
0.000 


0.978 
0.819 
0.235 
0.051 


0.987 
0.954 
0.618 
0.188 



Note. — Mean values for Pdetection for 
various period ranges (column 1) as function 
of observing run length (Fig. [SJ. Assumed 
parameters are given in th e caption of Fig. 
\3\ For discussion, see £13.21 



the night up between two or more fields, thereby decreas- 
ing the number of hours spent on each one of them. We 
illustrate the effect of such strategies in Fig. [4] and Ta- 
ble 2] in which we assume basically the same parameters 
as for Figures [3] and [5] (see £ )3.4p for purposes of com- 



parison. The situation shown in the bottom right panel 
can obviously only be achieved at numerically high lat- 
itudes on Earth during the respective winter season, or 
from space, but serves as a comparison to the scenarios 
encountered in transit surveys conducted from moderate 
latitudes. As in i j3.2[ the choice of strategy depends on 
the range of periods that is probed. 

One expected and visible effect in Fig. |4] is the de- 
creasing depth of spikes in Pdetection with longer lengths 
of night as the diurnal cycle becomes less of a factor in 
transit detection. As evidenced in Figures [7| and [TUl the 
spikes eventually disappear altogether when observing 
becomes uninterrupted, as, e.g., from space. 

3.4. Observing Cadence 

Observing cadence is primarily dependent on telescope 
and detector characteristics as well as target brightness, 
with the goals that a w is minimized, the target remains 
in the linearity regime of the detector, and the exposure 
time is not so long as to smear out phase information on 
any detectable planetary transit. 

Similar to Sj372]and §3.3[ however, the choice of cadence 
can also be used as a observing strategy parameter to in- 
crease the number of monitored stars at the expense of a 
lower sampling rate per field (by moving back and forth 
between fields between exposures, for instance). This 
effect is simulated in Fig. [5l and values for Pdetection 
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Fig. 3. — The influence of length of observing run on the detection efficiency, shown for 15, 30, 60, and 90 nights (indicated in the 
respective panel). cr w is assumed to be 5 mmag. The (blue) solid line indicates the detection efficiency o> = 0, and the (red) dashed 
for o> = 2 mmag. Additional parameters are: SNRt; ires / io ; < j = 7.0, 8 hours of observing every night, an observing cadence of 5 minutes, 
Rsta r = R@, Rplanet = 0.1-Rq, M sta r = Mq, M p i anet <g M@ . See Table[3]for mean values of Pdetection over various period ranges, and 
flOl for details. 



TABLE 4 

MEAN Pdetection VALUES FOR DIFFERENT 

Night Lengths 



Hours/Night 


3 


5 


10 


16 








= 0) 






0-5 days 
5-10 days 
10-20 days 
20-30 days 


0.949 
0.764 
0.358 
0.162 


0.970 
0.886 
0.543 
0.253 


0.989 
0.962 
0.846 
0.541 


0.997 
0.992 
0.979 
0.865 




(ay = 


2mmag) 






0-5 days 
5-10 days 
10-20 days 
20-30 days 


0.839 
0.182 
0.049 
0.017 


0.954 
0.462 
0.092 
0.028 


0.986 
0.912 
0.402 
0.075 


0.996 
0.985 
0.801 
0.237 



Note. — Mean values for Pdetection f° r var- 
ious period ranges (column 1) as function of 
length of night (Fig. |4j. Assumed parameters 
are given in t he caption of Fig. [4] For discus- 
sion, see i|3.3l 



for different period ranges are given in Table [5l which 
shows that the effect of changing from a cadence of one 
to several minutes does not greatly affect the calculated 
detection probability, especially for very small values of 



ay. It may therefore be worth considering changing be- 
tween fields every one or few exposures to increase target 
number. The effects of red noise produced for such an ob- 
serving strategy, however, such as flatfielding errors due 
to the fact that the stars may not be located in exactly 
the same position in the field as before, are dependent 
on aspects such as the pointing stability of the telescope 
used and would need to be explored for the respective 
observing setup. 

3.5. Minimum Number of Sampled Transits 
3.5.1. Two or more Transits 

Transit duration is dependen t on period (equation 3 in 
iSeager fc Mallen-Ornelasll2003l ). Observational cadence 
determines the number of expected data points in a single 
transit. The combination of number of nights in the ob- 
serving run and length of a given night sets the expected 
number of transits sampled during the monitoring cam- 
paign. 

As explained in §2.3[ we require a minimum of two 
transits sampled to constitute a detection. Note, how- 
ever, that other predictions of survey yields in the liter- 
ature use different numbers for different reasons, e.g., to 
be able to constrain period, which, with only two tran- 
sits detected, would be subject to significant aliasing un- 
certainties, depending on the time elapsed between the 
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Fig. 4. — The influence of the length of night on the detection efficiency, shown for 3, 5, 10, and 16 hours (indicated in the respective 
panel). a w is assumed to be 5 mmag. The (blue) solid line indicates the detection efficiency o> = 0, and the (red) dashed for oy = 2 mmag. 
Additional parameters are: SNRt^reshoid = 7.0, a 5-minute observing cadence, 60 consecutive nights, R s ta r = Rq, Rpianet = 0.1-Rq, 
M a tar = Mq, M p i anet -C Mq. See Table [4] for mean values of Pdetection over various period ranges, and 93.31 for details. 



TABLE 5 

MEAN Pdetection VALUES FOR DIFFERENT 

Observing Cadences 



Minutes 


1 


5 


30 


90 






(cr, 


= 0) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.984 
0.947 
0.758 
0.431 


0.983 
0.945 
0.753 
0.425 


0.974 
0.883 
0.535 
0.264 


0.725 
0.057 
0.021 
0.001 




(ay = 


2mmag) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.980 
0.864 
0.293 
0.071 


0.978 
0.819 
0.235 
0.051 


0.936 
0.214 
0.037 
0.000 


0.457 
0.015 
0.000 
0.000 



Note. — Mean values for Pdetection for 
various period ranges (column 1) as function 
of observing cadence (Fig. [5J|. Assumed pa- 
rameters are given in th e caption of Fig. [5] 
For discussion, see 93.41 



two s ampled transit events (see §7.4 in Ivon Braun et all 

Hoop). 

Fig. [6] and Table [6] show how the detection probabil- 
ity varies as a function of different minimum number of 
sampled transits. We note that, in this work, only one 



observation taken during a transit is enough to count this 
transit as sampled, but the detection is still a function of 
the transit SNR, as explained in §SJ as well as the number 
of sampled transits. 

It is interesting to observe that the detection probabil- 
ity of ay = 2 mmag approaches the o~ r — case for higher 
minimum number of sampled transits, showing how the 
" white noise only" case thus becomes increasingly equiv- 
alent to the realistic case with red noise present (Fig. 
[6] and Table [6] produce identical results for ay = and 
a r = 2 mmag for more than five or more transits sam- 
pled). In the absence of knowledge of oy, requiring at 
least three or four detected transits in the data could 
therefore serve as a alternative for calculating a conser- 
vative estimate of survey yield. 

3.5.2. Detections Based on Single Transits 

We now examine the case where it is deemed possible 
to detect a transit based on a minimum of one sampled 
transit. For this scenario, the detection probability may 
exhibit a non-intuitive behavior for very short and very 
long periods. To illustrate these points, we assume 5 
"nights" of uninterrupted 24- hour (e.g., space-based or 
polar ground-based) observing (see FiglZj) . 

Short periods imply short transit durations. For a 
given cadence, there are few points per sampled tran- 
sit. At the same time, short periods imply many tran- 
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period[days] period[days] 



Fig. 5. — The influence of observing cadence on the detection efficiency, shown for 1, 5, 60, and 90 minutes (indicated in the respective 
panel). cr w is assumed to be 5 mmag. The (blue) solid line indicates the detection efficiency a r = 0, and the (red) dashed for o> = 2 mmag. 
Additional parameters are: SNRthreshold = 7.0, 8 hours of observing every night, 60 consecutive nights, R s ta r = Rq, Rpianet = 0.1-Rq, 
M a tar = Mq, M p i anet -C Mq. See Table [5] for mean values of Pdetection over various period ranges, and i|3,4l for details. 



TABLE 6 

MEAN Pdetection VALUES FOR DIFFERENT 

Minimum Number of Sampled Transits 



Min. Number 


2 


3 


5 


10 






(oy 


= 0) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.983 
0.945 
0.753 
0.425 


0.981 
0.878 
0.339 
0.066 


0.965 
0.313 
0.023 
0.000 


0.549 
0.008 
0.000 
0.000 




(oy = 


2mmag) 






0—5 days 
5—10 days 
10-20 days 
20-30 days 


0.978 
0.819 
0.235 
0.051 


0.978 
0.819 
0.235 
0.051 


0.965 
0.313 
0.023 
0.000 


0.549 
0.008 
0.000 
0.000 



Note. — Mean values for Pdetection f° r vari- 
ous period ranges (column 1) as function of ob- 
serving cadence (Fig. [6J. Assumed parameters 
are given i n the caption of Fig. \E\ For discussion, 

see %nrn 



sits sampled for a given observing run length. With in- 
creasing orbital period, the number of sampled transits 
decreases as the number of data points per transit in- 
creases, though not at the same rate. This effect can 
be seen for the short period range in Fig. [7j and it is 



most intuitively understood by considering the case for 
oy = 0. Since a w — transit depth = 10 mmag (Fig. 
[7]), a detection simply requires SNI& reshold — 49 data 
points observed during transit (Eq. Qj) , even if they are 
all located in a single transit. Pdetection goes to zero at 
a period of 2.5 days, which, for the stellar and planetary 
masses and radii given in the caption of Fig. [JJ implies 
a transit length of around 161 minutes. For a 7-minutc 
observing cadence, 23 data points can thus be collected 
during a single transit. For a 2.5 day period, one would 
expect to have two transits present in a 5-day observing 
run (24 hours of observing per "night"), but two transits 
would only contain 46 data points, and not the 49 re- 
quired for SNRtransit > S TV Rthreshold- Pdetection Fig. 

[7J therefore goes to zero at that point. 

For periods longer than 3 days, however, the expected 
number of data points per transit is 24.5, and thus, sam- 
pling two transits somewhere in one's data is sufficient 
for a detection, resulting in an increase of Pdetection with 
period as it approaches 3 days. 

An inverse effect can be observed at periods longer 
than the duration of the observing run and is again most 
easily explained by considering the oy = case. As ex- 
pected, Pdetection drops to zero for periods longer than 
the observing run, since, although it is possible to sample 
one transit during the observing run, the length of this 
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Fig. 6. — The influence of the minimum number of transits sampled on the detection efficiency. Shown is the requirement that (at least) 
2, 3, 5, and fO transits (indicated in the respective panel) be sampled. cr w is assumed to be 5 mmag. The (blue) solid line indicates the 
detection efficiency oy = 0, and the (red) dashed for o> = 2 mmag. Additional parameters are: SNRj^reshold = 7.0, a 5-minute observing 
cadence, 8 hours observing per night, 60 consecutive nights, R s tar = Rq, Rpianet = 0.1-Rq, M s tar = M@, M p i ane t <C Mq. See Table[6] 
for mean values of Pdetection over various period ranges, and ^13. 5.11 for details. 



one transit is too short for enough data points (49) to be 
sampled during its duration (and it is obviously impos- 
sible to sample more than one transit, as indicated by 
the dotted line). This situation, however, changes as the 
period approaches 24 days (transit duration = 343 min- 
utes). From that point on, the transit duration will be 
long enough to fit 49 data points at a 7-minute cadence. 
Therefore, if a single transit is observed, it is long enough 
to gather enough data to fulfill the threshold SNR crite- 
rion. Ultimately the probability that any transit occurs 
at all during the observing run, and thus the detection 
efficiency, approaches zero. 

3.6. Transit Depth and Duration 

Throughout the paper, we calculate the tran- 
sit depth and duration accor ding to equations in 
ISeager fc Mallen-Ornelai 41)03). thereby assuming a 
central transit, i.e., i = 90° and b = 0, as well as so- 
lar and Jupiter values for stellar and planetary radii and 
masses. In this subsection, we set those parameters to 
certain values that may not be consistent with physical 
laws, but are meant to illustrate the behavior of the de- 
tection probability as a function of transit duration and 
depth. 

Transit depth is primarily a function of stellar and 
planetary radius, which we set to solar and Jupiter val- 



ues before, resulting in a transit depth of 0.01 mag or 
1% of the relative flux. In order to show how much more 
challenging the detection of smaller planets is, or con- 
versely, how much easier the detection of larger planets 
is for given observing parameters, we vary transit depth 
in Fig. [5] and show mean Pdetection values in Table [7J 

For a parent star with R = Rq, the panels repre- 
sent planets of 0.3 Rj up iter (top left), 0.7 Rj up iter (top 
right), 1.0 Rj up iter (bottom left), and 1.4 Rj up iter (bot- 
tom right). Note that an Earth-sized planet would have 
a radius of around 0.1 Rjupiter and produce an eclipse 
around a solar-sized star with depth of 10% of that as- 
sumed in the top left panel, i.e., 0.0001 mag. It is worth 
pointing out how significant the difference in detection 
probability is for shallow transits between ay = and 
a r 7^ 0, substantiating the claim that space-based ob- 
serving is necessary to find very small planets (see also 
£ j4.1[) , as recently evidenced by the discovery of CoRoT- 
7b (Leger et al. in preparation; Rouan et al. in prepara- 
tion; Bouchy et al. in preparation). 

Transit duration is a function of orbital period, i, and 
stellar and planetary masses and radii. Rather than fol- 
lowing the physical dependence on period, we set transit 
duration to fixed values of 1, 2, 5, and 10 hours in the four 
panels of Fig. [9] (see also Table [8]), thereby still assuming 
values for the parameters mentioned in the figure cap- 
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Fig. 7. — The behavior of the detection probability as a function of period for a short, white-noise dominated monitoring campaign, tj w 
is assumed to be 10 mmag. oy = for the solid (blue) and dotted (black) lines; oy = 0.5 mmag for the dashed (red) line. The solid (blue) 
and dashed (red) lines no longer assume a minimum number of sampled transits, whereas the dotted (black) line requires a minimum of 
two transits to be present in the data. Additional parameters are: Subthreshold = 7.0, a 7-minute observing cadence, continuous observing 
(i.e., 24 hours observing per "night"), 5 consecutive nights, Rstar = Rq, Rpianet = 0.1- Rq, M s tar = Mq, M p i anet = Mj upiter . The inset 
is a zoomed display of the behavior of P, detection f° r periods between 1 and 5 days. See i|3. 5.21 for discussion. 



TABLE 7 

Mean Pdetection Values for Different Transit Depths 



Transit Depth 


1 mmag 


5 mmag 


10 mmag 


20 mmag 




(ov = 0) 


0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.009 
0.000 
0.000 
0.000 


0.977 
0.916 
0.677 
0.355 


0.983 
0.945 
0.753 
0.425 


0.984 
0.947 
0.755 
0.428 




(o"j 


= 2mmag) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.000 
0.000 
0.000 
0.000 


0.435 
0.005 
0.000 
0.000 


0.978 
0.819 
0.235 
0.051 


0.983 
0.947 
0.755 
0.428 



Note. — Mean values for Pdetection f° r various period 
ranges (column 1) as function of transit depth (Fig. 0. As- 
sumed paramet ers are given in the caption of Fig. [8| For 
discussion, see i|3.6l 



tion, including a transit depth of 0.01 mag. While longer 
transits are obviously easier to detect, the increase of de- 
tectability with transit duration is slow but sensitively 
dependent on a r . 



TABLE 8 

MEAN Pdetection VALUES FOR DIFFERENT TRANSIT 

Durations 



Transit Duration 


1 hour 


2 hours 


5 hours 


10 hours 






(oy 


= 0) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.976 
0.896 
0.506 
0.193 


0.981 
0.920 
0.585 
0.240 


0.989 
0.957 
0.776 
0.409 


0.996 
0.987 
0.938 
0.690 




(oy = 


2mmag) 






0—5 days 
5-10 days 
10-20 days 
20-30 days 


0.964 
0.414 
0.058 
0.010 


0.975 
0.676 
0.126 
0.028 


0.986 
0.879 
0.272 
0.050 


0.995 
0.974 
0.609 
0.121 



Note. — Mean values for Pdetection f° r various period 
ranges (column 1) as function of transit duration (Fig. [9t . 
Assumed para meter s are given in the caption of Fig. [5] For 
discussion, see £13.61 



4. APPLICATION AND EXAMPLES 

The examples used in Sj3] to illustrate the influences 
of various observing strategy and astrophysical param- 
eters on the transit detection probability resemble ob- 
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Fig. 8. — The influence of the depth of transit on the detection efficiency, shown for 1, 5, 10, and 20 mmag (indicated in the respective 
panel). a w is assumed to be 5 mmag. The (blue) solid line indicates the detection efficiency o> = 0, and the (red) dashed for oy = 2 mmag. 
Additional parameters are: SNRt^reshoid = 7.0, a 5-minute observing cadence, 8 hours observing per night, 60 consecutive nights, transit 
duration calculated using solar radius and ma ss fo r the star and Jupiter values for planetary mass and radius. See Tabic 171 for mean values 
of Pdetection over various period ranges, and A3, 61 for details. 



serving campaigns typical of the very successful wide- 
field transit surveys such as HAT, TrES, XO, SWASP 
etc (e.g., iBakos et all I2007t 10 'Donovan et all 120061 : 
iMcCullough et alj 120061 : iPollacco et alj 120061 ). In con- 
trast, this section shows examples and consequences of 
observational window functions for fundamentally differ- 
ent setups of monitoring projects. 

4.1. Space Based Surveys 

Compared to ground-based c ounterparts, space- based 
transit surveys such as CoRoT (Ba glin et alJ l2006) have 
the principal two advantages that (a) they are not subject 
to interruptions in observing due to the diurnal cycle 
(see §3.3p . and that (b) they do not need to deal with 
the Earth's atmosphere (see fl3.1|) . The latter aspect in 
particular makes them the currently only realistic option 
of detecting Earth-sized planets around sun-like stars, 
which is one of t he explicit goals of t he recently launched 
Kepler Mission (jBorucki et al.ll2009D . 

Fig. [10] shows the detection probability for simulated 
space- based surveys of various lengths, loosely modeled 
after the long and short observing runs by the CoRoT 
satellite, thereby assuming somewhat generic parameters 
for survey strategy and photometric precision (see cap- 
tion). The solid and dashed lines respectively indicate 
the detection probabilities for a Jupiter-sized planet and 



for the recently discovered exoplanet CoRoT-7b around 
its parent star, a KO dwarf (transit depth ~ 0.5 mmag; 
period ~ 0.9 days). We note that our simulations of an 
Earth-sized planet around a solar-type star produce a 
detection probability of zero for all periods. 

4.2. Synoptic Surveys 

Synoptic surveys typically provide high-quality pho- 
tometric time-series data of very low cadence but over 
extended periods of time. Thus, they are not pri- 
marily designed to find planetary transits but never- 
theless present data sets that are worth prob i ng fo r 
their existence (see for instance iPlavchan et alJ [2008) . 
In fact, several transiting planets have been discov- 
ered a po steriori in the Hipparcos a rchives such as HD 
209458b dRobichon fc Arenoul |2000D an d HD 189733b 
(jHebrard fc Lecave licr Dcs Etangsll2006[ ). 

The panels in Fig. [TT] are produced by observational 
window functions of synoptic surveys loosely based on 
the (the future, ground-based) La rge Synoptic Survey 
Telescope fLSST: llvezic et aT1l2008l ) in the top panel and 
(the space-based) Hipparcos mission (|Perrvman fc ESAl 
Il997h in the bottom panel. For both panels, we require 
a S N Rthreshoid = 7.0, at least two sampled transits, and 
assume solar and Jupiter values for stellar and planetary 
mass and radius. The principal differences between the 
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Fig. 9. — The influence of transit duration on the detection efficiency, shown for 1, 2, 5, and 10 hours (indicated in the respective panel). 
cr w is assumed to be 5 mmag. The (blue) solid line indicates the detection efficiency a r = 0, and the (red) dashed for o> = 2 mmag. 
Additional parameters are: SNR^reshoid = 7.0, a 5-minute observing cadence, 8 hours obser ving per night, 60 consecutive nights, transit 
depth = 0.01 mag. See Table|8]for mean values of Pde.tec.tion over various period ranges, and A3, 61 for details. 



two window functions in Fig. [TT] are due to the different 
assumptions in a w (5 mmag for top panel; 1.5 mmag 
for bottom panel) and oy (1 mmag for top panel; 0.5 
mmag for bottom panel), and the different number of 
data points obtained over different lengths of time. 

For the top panel, we assumed a cadence of a single, 
30-second exposure time image every three nights, ac- 
cumulated over around eight years, such that the total 
number of images is 1,000. The mean value for Pdetection 
over various period ranges are as follows: 

• 0-10 days: < Pdetection > = 0.555; 

• 10-50 days: < Pdetection > = 0.239; 

• 50-100 days: < Pdetection > = 0.025; 

• 100-200 days: < Pdetection > = 0.007; 

For the bottom panel, we chose a cadence based on 
the actual observations of a Hipparcos star with 190 
epochs, downloaded from the NASA Star and Exoplanet 
Database 2 . Basically, the 190 observations were obtained 
over three years in groups of several images every few 
tens of days. The mean value for Pdetection over various 
period ranges are as follows: 

2 http://nsted.ipac.caltech.edu 



• 0-10 days: < Pdetection > = 0.746; 

• 10-50 days: < Pdetection > = 0.156; 

• 50-100 days: < Pdetection > = 0.034; 

• 100-200 days: < P d etecUon > = 0.012; 

Finally, it should be noted that we ran the equivalent 
simulations to the ones in Fig. [Til but thereby assum- 
ing an Earth-sized planet instead of a Jupiter-sized one. 
Both detection probabilities were identical to zero for all 
periods. 

5. DISCUSSION AND CONCLUSION 

This work quantitatively illustrates the influence of a 
number of observing strategy and astrophysical parame- 
ters on the detection efficiency of existing planetary tran- 
sits as a function of period under general assumptions 
listed in ^2.11 and parameters given in the various figure 
captions. The influences of red and white noises upon 
this detection efficiency are first examined in their own 
right, and then included in every simulation of the afore- 
mentioned parameters. Red noise is confirmed to be the 
dominant challenge to overcome in the search for plane- 
tary transits, as seen in the discussion in H3.ll Figures [JJ 
and[2] and Tables [JJ and [2j All parameters being equal, a 



Observational Window Functions in Planet Transit Surveys 



13 



IX 



0.8 



0.6 



0.4 



0.2 



£ 1 



0.8 



0.6 



0.4 



0.2 



1 

0.8 
| 0.6 

V 

o7 0.4 




M 1 1 1 1 I 1 1 1 1 limp 



0.5 1 1.5 2 2.5 3 
period[days] 



I I I =F 



150 



2i)0 



100 
period[days] 



150 



200 



Fig. 10. — Simulated space survey loosely based on the CoRoT short runs (bottom panel) and long runs (top panel). The solid (red) 
line indicates the detection probability of a Jupiter-sized planet around a solar-type star, the dashed (blue) line indicates the detection 
probability of the recently discovered exoplanet CoRoT-7b around its (K0 dwarf) parent star (transit depth ~ 0.5 mmag; period ~ 0.9 
days). The inset in the top panel represents a zoomed display of the 0.5-3.0 day period range. Parameter values assumed in this simulation 
are: SNR t f lres i lo i ( i = 7.0, o~ w = 1 mmag, a r = 0.5 mmag, continuous observing with a 15-min cadence over 30 days (bottom panel) and 
150 days (top panel). See £|4.1l for details. 



factor 4 increase in ay produces a much more significant 
reduction in Pdetection than a factor 10 increase in a w . 

In particular, we explicitly address controllable strat- 
egy parameters such as the number of nights for which 
one may choose to monitor a given target held, the num- 
ber of hours per night one may stay on this field, and the 
observational cadence with which the field is monitored. 
We furthermore examine the influence of astrophysical 
parameters on detection efficiency, such as transit depth 
and duration. Finally, we look at parameters typically in- 
volved in the calculation of the projected yield of a given 
transit survey such as the minimum number of transits 
required for detections, and illustrate two non-intuitive 
effects that occur when the criterion of a minimum num- 
ber of sampled transits is abandoned and detection is 
based only on SNR. Along with visualization of the ef- 



fects caused by the various parameters in the figures, we 
provide quantitative means of comparison for different 
period ranges in the accompanying tables. 

A consideration that did not factor into the calcula- 
tion of Pdetection is the fraction of data points outside 
of transit (we assumed that this fraction is much higher 
than the number of points sampled in transit). In order 
to detect a transit in one's data, one needs to have both 
brightness levels well measured such that the difference 
between them becomes signficant enough to enable a de- 
tection. For instance, an observing run that only obtains 
data during transit would, by the metrics used in this pa- 
per, detect the transit, provided the SNR is high enough. 
In real life, however, the data would appear perfectly 
flat and no sensible algorithm would flag the signal as 
a possible planetary transit. Obviously, sparse cadences 



14 



von Braun et al. 




period[days] 



Fig. 11. — Simulated synoptic surveys loosely modeled after LSST (top panel) and Hipparcos (bottom panel). The insets are zoomed 
regions for the shorter periods. We assume solar-type stars and Jovian planets in these simulations. The 1,000 obse rvat ions in the top 
panel are obtained over eight years, whereas the 190 observations in the bottom panel are taken over three years. See £|4.2l for details. 



are more susceptible to this admittedly pathologic pitfall 
than well sampled ones. One scenario where one may 
encounter a problem like this would be in the attempt to 
detect transits among long-period planets discovered by 
radial velocity work, which would tend to exhibit long 
transit durations. 

More generally, we caution that the detection of an ex- 
isting transit in "real life" is dependent on a large number 
of properties of the data reduction and analysis pipeline 
and transit detection methods, including human expe- 
rience and human error potential, which cannot possi- 
bly be parametrized as an ensemble or included in any 
code. Therefore, the significance of our results and pre- 
dictions, although quantitative, are necessarily subject to 
an unknown fud g e or s caling factor. As pointed out by 
iBeattv fc Gaudl ([2008), the non-uniformity of the defi- 
nitions of detection criteria cause the largest uncertainty 
in transit survey yield predictions. 

Nevertheless, we specifically allowed for parameters 



that are typically calculable in transit survey designs 
to be used as input to the code in order to make it as 
practically applicable as possible. Consequently, even 
in the presence of the unknown fudge factor mentioned 
above, comparisons between different observing strate- 
gies is quantitatively possible to optimize survey yield. 
For instance, under some circumstances it appears much 
more favorable to increase the observational cadence 
to add a second monitoring field to one's project than 
switching fields in the middle of the night or in the middle 
of the observing run (see §j j3.2[l3~3l and l3.4j) . provided it 
is possible to repeatedly achieve very good pointing of the 
telescope to reduce the additional red noise component 
that might arise from flatfielding errors otherwise. Thus, 
it may well be advisable for transit surveys to consider 
this trade-off between cadence and fields monitored as it 
can lead to a dramatic change in the predicted planet 
yield of the survey. 
Furthermore, the examination of the effects of red and 
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white noise in £)3.1l and throughout the paper give quanti- 
tative insight into what size of planet one may expect to 
realistically detect in one's data, given observing strat- 
egy parameters. In general, the depth of transit one may 
hope to detect in one's data needs to be larger than the 
magnitude of ay, as seen in §3.61 and evidenced Figures [8] 
and [TU] This is confirmed very well in Fig. [IJjJ showing 
that the detection of CoRoT-7b around its parent star, 
given their sizes and orbital period, are right at the limits 
of the CoRoT satellite for a single long run (i.e., without 
combining data from several runs). 

The code used for all calculations in this paper is avail- 
able from KvB upon request. 
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